Mapping the Dark Energy with Varying Alpha 
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Cosmological dark energy is a natural source of variation of the fine structure constant. Using 
a model-independent approach we show that once general assumptions about the alpha-varying 
interactions are made, astronomical probes of its variation constrain the dark energy equation of 
state today to satisfy —1 < Wf < —0.96 at 3-sigma and significantly disfavour late-time changes in 
the equation of state. We show how dark-energy-induced spatial perturbations of alpha are linked 
to violations of the Equivalence Principle and are thus negligible at low-redshift, in stark contrast 
to the BSBM theories. This provides a new test of dark energy as the source of alpha variation. 



The universe appears to be accelerating in response 
to a dark energy whose identity remains a mystery. We 
show how the current astronomical data supporting a 
5a detection of time variation in the fine-structure 'con- 
stant', a — e 2 /he, at redshifts 0.2 < z < 3.7 [1,2] can pro- 
vide direct evidence distinguishing a dynamical dark en- 
ergy from a cosmological constant, which is very challeng- 
ing for standard cosmological tests, e.g. [3]. Unlike cos- 
mological tests which depend on the integrated properties 



of the dark energy, measurement of 



Aq(z) _ a(z)-a(0) 

^oy— 

would allow us to map the dark energy as a function of 
z, see also refs. [4-6]. 

Mapping the dark energy - We will consider a varia- 
tion of a generated by a non-trivial gauge kinetic function 
(GKF), Zp((f)), leading to an electromagnetic lagrangian 
jZp(<j})F, J , u F IJ ' v , where <j> is some scalar field and F^ v the 
usual Maxwell tensor. (We expect Zp ^ 1 in string the- 
ory through the coupling to the dilaton [7] or through 
non-renormalisable interactions in supergravity [8]) In 
this case a(z) cx Zp(4>)~ 1 so if we know Zp and Aa(z)/a 
precisely we can deduce <fi(z). 

Is it natural for <f> to be the dark energy? The inter- 
action Zp((f))F 2 is not renormalisable in 4d and hence 
Zp = Zp(cj)/M) where M is the mass scale at which 
the effective theory description (usually supergravity) 
breaks down. If we assume M is the Planck mass, 
Alpi ~ 10 19 GeU, only very large expectation values of 
<t> will lead to measurable effects. This is what hap- 
pens in many dark-energy models that typically require 
large field values, <p/M p i = 0(1) today to satisfy the 
slow-roll condition, M^V^ < V, on the potential V(4>) 
needed for an accelerating universe. We will show, using 
a model-independent approach, how one can constrain 
the dark-energy's equation of state. Note that such mod- 
els for varying a differ from those of Bekenstein-Sandvik- 
Barrow-Magueijo (BSBM) [9] in which a-variations are 
driven by coupling to the charged non-relativistic mat- 
ter alone and so the associated scalar field cannot be the 
dark energy. 

We can also consider spatial variation of a [10,2]. Dark 
energy should cluster on scales larger than its Compton 



— 1/2 

wavelength A c oc , which exceeds 100-Mpc in stan- 

dard quintessence models [11], although in more exotic 
models this need not be true [12]. A map of the spatial 
variation of a can potentially probe both the power spec- 
trum and Compton wavelength of fluctuations in the dark 
energy. Consider the spatial fluctuation 5a = a(x^) — a. 
Quantities with an overbar are evaluated in the back- 
ground cosmology at the same event. To first order in 5(f) 
we have 



<$a(x,i) 



(d In a M 



(1) 



Evaluated today, this is directly related to Weak Equiva- 
lence Principle (WEP) constraints [13] on varying-a the- 
ories which, assuming large A c so that we may use the 
background solution for cj) even in the nonlinear regime* , 
imply 5a/a\o < I0~ 5 5<fi/M p i, today. Since we expect 
5(f> <C M p i, spatial fluctuations of a should be negligible 
at low z if dark energy (or any light scalar field) is the 
source for the time-variation of a. This is in agreement 
with current data [2]. But at high z the spatial fluctua- 
tions of a could be significant. Eq. (1) implies that the 
power spectrum of a fluctuations, Pg a (k) = ^p-|<5ckfc| 2 is 
anti-correlated with the spectrum of dark-energy fluctu- 
ations, Ps(f,(k). Their ratio, which we call the a-bias b a , 
will in general be time-dependent with potentially inter- 
esting implications for the CMB [16]. 

Zp degeneracy - Our discussion has been predicated 
on the (unreasonable) assumption that we know the func- 
tion Zp. A standard approach [5] is to assume slow roll 
for 4> in some potential V(<p) and expand Zp around 4>q, 
the value of cf> today (see for example ref. [6]). The leading 
term gives Aa(z)/a oc A(f>/M p i. However this is overly 
restrictive since Zp may not be analytic in <j), the dark en- 
ergy may not derive from a slowly-rolling scalar field [15], 



*By contrast, in BSBM varying-a theories spatial varia- 
tion in a is strong because it is driven by gradients in 
electromagnetically-charged matter, removing the power of 
local constraints [14]. 
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(a) Two different GKFs 



(b) Degeneracy of ty and z, 
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„=-!, z=0.6, $,,=0.31 
,=-0.985, z=5, <t>„=0.01 



FIG. 1. Degeneracies: (a) Aa(z)/a for the same 
dark-energy model with two different choices of Zp. Note 
that Zf = A sin cj> allows both positive and negative values 
of Aa(z)/a. (b) The z t — <fio/M degeneracy in fitting the 
Aa(z)/a data. Both models have almost the same y/ 2 value. 
The WEP data breaks this degeneracy. 



or Zp may change very rapidly with cf>, e.g. Zp oc so 
even if <fi ~ ln(i), Aa(z)/a would exhibit power-law vari- 
ation with redshift. Hence there is a perfect degeneracy 
between Zp and the dark-energy dynamics, shown in Fig 
(la) by plotting Aa(z) / a for the same dark-energy model 
with two very different GKFs: Zp((f>) = \(f> and Asin(^). 
If the resulting curves were assumed to arise from the 
same GKF, they would have been interpreted as coming 
from two completely different dark-energy models. 

To further illustrate this degeneracy, consider the left 
branch of the Oklo analysis of the neutron capture reso- 
nance [19] which suggests Aa/a = (0.88 ± 0.07) x 10~ 7 . 
This corresponds to a larger a at z = 0.13, in contrast 
to the Webb et al data which strongly suggests a was 
smaller at z > 0.5 with Aa/a = -(0.543±0.116) x 10~ 5 , 

It has been claimed that this is inconsistent with a 
slowly-rolling quintessence origin; true if Zp is monotonic 
in cj). However, a GKF such as Zp((f>) = Asin(0) can give 



both positive and negative 



Aa(z) 



and match the appar- 



ently inconsistent data with monotonic (j)( z )> as shown 
in Fig. 1(a). In fact, given any data Aa(z)/a and any 
smooth 4>(z) we can reconstruct a corresponding GKF 
which will exactly fit the data by defining the GKF to be 
Z F {(j){z)) = a*/a(z) - 1 where a* ~ (1 + Z F (0 O ))/137 
normalises a(z). 

Dark energy can induce variation of a in other ways. 



We can choose a tensor, 



■Jem 



to raise and lower indices 



on that differs from the spacetime metric g^ v used 
to form the Ricci scalar R [17]. A simple example is 



9em = ST 



(V0)" 
M 2n+2 



(2) 



where V M is the covariant derivative w.r.t. g^ v . This 
yields Z F = Z F (\ V^|/A/ 2 ), so Aa/a ~ {H/M). If 
M = M p i such interactions are negligible today; but 
if M ~ 10~ 3 eV, the models can easily fit the data at 
z < 3.5; however, if d> is a tracker field then d> — > oo 



as z — > oo. If Zp{4>) oc 1 + (<j)/M 2 ) n for any n, then 
\Zp(oo)\ —> oo and a(z) — > rapidly. In fact, we have 
found that models of this form that match the data at 
z < 3 but strongly violate the nucleosynthesis constraint 
|^(10 10 )| < 0.1 [18]. This can be avoided by choosing 
a GKF with Zp{<f> — > oo) = constant, or a model with 
bounded |(/>|, as in condensation models of dark energy 
[12], but this goes beyond the scope of the paper. t. This 
shows the importance of theoretical input into choosing 
Zp and the need for general parametrisations of Zp. Be- 
low we investigate to what extent we can expect to dis- 
tinguish between different Zp and different dark energy 
equations of state. 

Consistency - Before we discuss mapping of the dark 
energy we address in more detail the suggestions in the 
literature (e.g. [4,20]) that it is difficult, or even impos- 
sible, for a slowly-rolling scalar field to match both null 
geonuclear data [19,25] at z < 0.5 and the quasar data at 
z > 0.5 in Zp theories (it is easier in BSBM theories be- 
cause there is strong spatial variation [14]). To quantify 
the extent to which dark energy can match the data we 
use model-independent approach [12,21-23] which uses 
w(z) = p$/ p$ instead of choosing a particular potential. 
An expedient choice is 



w(z) =vjf[l + exp(((z - z t )/A)\ 



(4) 



so that w(z — > oo) = 0; Wf is constant; z t is the redshift 
at which a transition to dark-energy domination occurs, 
while the constant A controls its width. This form accu- 
rately describes most quintessence models and some ex- 
otic alternatives such as vacuum metamorphosis [15]. We 
fix A = zt/30 as in [21,12] so w(0) = Wf, and perform 
a likelihood analysis over the free dark-energy parame- 
ters (wf,z t ) and the parameters of the GKF. We have 
already seen that if Zp is completely free one can match 
any data with any dark-energy dynamics. Therefore, we 
need theoretical input or a general parametrisation of the 
GKF. In the spirit of the model-independent approach we 
consider an expansion of the GKF in the form: 



Zp = 1- T n ip n , r„ 



Mb 



(5) 



where ip(z) = <f>(z)/<po and r„ > are dimensionless 
couplings. We limit ourselves to n < 4 for computational 
reasons and consider each n separately. This choice of 
Zp((j>) expands our set of parameters to (wf,Zt',(3,jj,ri). 



^Another example, the 1-loop QED lagrangian in curved 
spacetime, [24] yields a GKF of the form 



Z F = 1 + aiT + bR^R^ + cR^sR 



(3) 



where a, b, c are constants and i?^ 7 * is the Riemann tensor. 
However, these terms are Planck suppressed and are not ex- 
pected to be influential today. 
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FIG. 2. The individual best-fits for each value of n. The 
difference in Aa/a between the models for z < 3 is less than 
2 x 1CF 6 and yet Figure (3) shows that the combined data 
strongly favour n — 3. This reflects the highly non-Gaussian 
nature of the x 2 hyp ersurf aces. 



The best-fits for each n, varying the remaining parame- 
ters, are shown in Fig. (2). 

Our first Aa(z)/a data set is {QSO}, the quasar 
absorption data [1,2] providing the evidence for time- 
variation of a. The second set of data is {QSO + Oklo 
+ Re}. {Oklo} is the null branch of the Oklo data 
(|^f (z = 0.13)| < lCT 7 ) [19]; {Re} is the from the Os/Re 
ratio from radioactive decay 187 Re — > 187 Os in iron mete- 
orites, [25] (l^f = 0.44)| < 3xl0" 7 ) close to the lowest 
redshift quasars. Last, we add composition-dependent 
WEP constraints on the difference in the acceleration, 
a, of different bodies A and B, rj = (Aa/a)^ < 10~ 12 
[13,4] that arises because of their different numbers of 
charged nucleons. Allowing for gauge coupling unifica- 
tion, r\ oc (d\na/d(j>) 2 [7] (and note eq. (1)) and this 
translates into the bound 



(i - r. 



M 



< io- 



Mr,! 



(6) 



today. 

The best-fit reduced-x 2 for the three data sets are 
4.7/8, 4.8/10 and 10.4/11 respectively. Adding the WEP 
constraint clearly has the largest effect but the \ 2 Y3 ^~ 
ues show that general dark-energy models are compatible 
with the current a{z) data. 

Closing the circle - Is it possible to gain information 
about Zp(4>) and the dark energy simultaneously? To 
answer this question we compute the 1-d likelihoods for 
each one of the parameters in the set (wf,Zt\/3,j%,ri) by 
marginalising over the remaining parameters on a grid 
consisting of 2.3 million points using the following precise 
expression for Aa(z)/a: 



Aa(z) 1 -¥>"(*) 

— J- n" 



a 'l-r„¥>»(2f) 
Our results, are shown in Fig. (3). 
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FIG. 3. The marginalised ld-likelihoods for n, f3, <f)o/M and 
zt for each of the three combinations of data sets. The curves 
represent the final likelihoods including all the data. Note 
how the data strongly favour n — 3 and how the WEP data 
selects 4>o/M ~ 10 -2 as the most favoured value. 



In particular, a dimensionless coupling (3 ~ 10 -3 — 10 -5 
is preferred by all the data. The role of n is interesting. 
While for each n it is possible to find sets of parameters 
with nearly equal x 2 > Fig- (2), the marginalised like- 
lihood for n is selective, Fig. (3). This reflects the fact 
that the likelihood hypersurface is highly non-Gaussian, 
with a large number of nearly degenerate minima. How- 
ever the most likely values of our parameters depend 
on the data set used, casting some doubt on the ma- 
turity of current data. Considering only the {QSO}data 
set, n = 1 is preferred. However, when the low-redshift 
{Oklo/Re} and WEP constraints are added then n = 3 
is significantly favoured, which is a promising sign for fu- 
ture determination of the GKF. If we consider only the 
{QSO + Oklo + Re}data, there is an interesting degener- 
acy between dark-energy models with low (high) z t and 
high (low) values of 4>o/M. This is illustrated in Fig. 
(lb) by two good fits to the data with parameter values 
(z t ,(j) /M) = (5.0, 10" 2 ) and (0.6, 10 05 ) respectively. 
In order to map the dark energy we calculate 4>{z) = 
— Jq 4>(z)dz/(H(l+z)) where 2 (z) = (l+w(z))p c f > and 
P4> (z) = p exp(3 f(l+w)dz/(l+z)) [21,12]. We set fi = 
0.7 today. These equations imply that as uif — > — 1 and 
Aa/a — ► for z < z t . Clearly, the tighter the null results 
on Aa at low redshift become, the closer wf is pushed 
to —1. This can be clearly seen in the likelihood curve 
for Wf in Fig. (4) when the {Oklo +Re} and WEP data 
are added. The WEP constraint is particularly useful 
since it breaks the degeneracy between (3, <pQ and M in 
eq. (7) which only depends on the combination T n . The 
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constraint is weakest if we choose M = M v \, as we do. 
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FIG. 4. The 1-d likelihood for log(l + w f ). The dashed 
line corresponds to the quasar data, the dot-dashed line to 
the quasar+Oklo+Re-Os data, and the solid line to all the 
data (including the WEP violation constraint). The full data 
set gives —1 < Wf < —0.96 at 3<j. 

The complete set of data favour both wf close to — 1 
and z t > 2.5; see Figs (4) & (3). {Oklo + Re} favours 
4>(z < 0.45) = and hence Wf ~ — 1. Adding WEP 
restricts —1 < Wf < —0.96 at 3er with wt = —1 signifi- 
cantly disfavoured, see Fig. (4). This is because Wf = — 1 
would imply Aa(z)/a = for z < z t , and to match the 
quasar data would then require z t < 0.5, inconsistent 
with the likelihood curve for z t . Small values of z t values 
are disfavoured because, as shown in Fig. (1), there is 
a z t — (<po/M) degeneracy in matching the quasar data. 
Models with small z t require large values of (f>o/M, but 
large <po/M is disfavoured by the Oklo/Re data (Fig 3). 
Conversely, very small 4>o/M are disfavoured (at constant 
r n ) by the WEP data, eq. (6). 

For n = 1 (the standard dimension-5 case [4,5]) the 
constraint (6) yields ft < 10~ 5 (M/M p i), while for n > 2 
it is a joint constraint on T n and <f> /M. In the next 
decade the STEP [26], GG [27] and MICROSCOPE [28] 
experiments promise sensitivities up to t] < 10 -18 and 
hence a null detection would imply d\na/d(j) < 10~ 8 to- 
day, virtually ruling out n = 1. To investigate potential 
future constraints we can impose this expected bound on 
our best-fits for each n. The resulting reduced-x 2 s all 
exceed 38/11 and would all be excluded by this future 
data. 

In summary, by using a model-independent approach 
to the dark energy, we have shown that the current con- 
straints on Aa(z)/a can be well-matched by dark-energy 
models. For a large family of varying-a theories, we 
have shown how spatial fluctuations in a are correlated 
with fluctuations of the dark energy but need to be neg- 
ligibly small at low-redshifts due to Weak Equivalence 
Principle constraints, so providing a clear test of the 
scenario. We have pointed out that the gauge kinetic 
function, Zp, responsible for the variation of a is com- 
pletely degenerate with the dark energy dynamics so that 
any dynamical dark-energy model can be made to fit the 
Aa ^ data. Nevertheless, with a reasonable parametri- 
sation of Zp the current data already yield interesting 



constraints on Zp. When the WEP constraint is added 
the current equation of state is forced to lie in the range 
— 1 < Wf < —0.96 at 3(7, the tightest constraints on 
the dark energy yet found. This is unlikely to change 
much if one varies other cosmic parameters. Late-time 
change in w(z) at z < 2 and ^ > 10 _1 are also strongly 
disfavoured, in contrast to standard CMB/LSS results 
[21,12], raising the possibility that a complete combined 
analysis will rule out most dark-energy models or even 
perhaps the paradigm of dark-energy-induced variation 
of a. 

Finally, we note that the GKFs we have considered 
break the conformal invariance of the Maxwell equations 
and may thus allow the seeding of magnetic fields dur- 
ing inflation. This would provide an interesting way of 
linking the early universe to observations in the present 
epoch. 
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